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ABSTRACT:  The  history  of  magnetic  circular  dichroism  (MCD)  spectroscopy  in  the 
study  of  alkali  metal /rare  gas  (M/Rg)  cryogenic  systems  is  reviewed  in  the  context  of 
developing  a  better  understanding  of  alkali  metal/hydrogen  systems  of  current  interest 
to  the  U.  S.  Air  Force  as  enhanced-performance  cryogenic  rocket  propellants.  A  new 
theory  for  simulating  the  MCD  spectra  of  M/Rg  systems  is  presented  together  with  a 
careful  discussion  of  the  theory's  implicit  and  explicit  approximations  and  their 
implications.  This  theory  uses  a  classical  Monte  Carlo  (MC)  simulation  scheme  to  model 
the  perturbing  effects  of  the  Rg  environment  on  the  2S  — *  2P  MCD-active  transition  of 
the  M  atom.  The  theory  sets  up  the  MC-MCD  simulation  as  a  6  X  6  matrix  eigenvalue/ 
eigenvector  problem  in  the  2P  manifold  in  which  are  included  the  effects  of  M-Rg 
interactions,  metal  atom  spin-orbit  coupling  in  the  2P  manifold,  magnetic  Zeeman 
perturbations  of  the  2S  and  2P  manifolds,  Boltzmann  temperature  factors,  and  electric 
dipole  transition  moment  integrals  for  left  circularly  polarized  (LCP)  and  right 
circularly  polarized  (RCP)  light.  The  theory  may  be  applied  to  any  type  of  trapping  site 
of  the  host  M  in  the  guest  Rg  matrix;  a  single  atom  substitutional  metal  atom  trapping 
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site  (one  host  Rg  atom  is  replaced  by  one  guest  M  atom)  is  modeled  in  this  study  for 
M  =  Na  and  Rg  =  Ar.  Two  temperature  factors  are  used  in  these  simulations;  a  lattice 
temperature  to  model  the  mobility  of  the  Rg  lattice  and  a  magnetic  temperature  to 
model  Boltzmann  factors  in  the  2S  ground  manifold.  The  6X6  eigenvalue /eigenvector 
problem  is  solved  for  a  number  of  randomly  generated  and  suitably  averaged  Rg 
configurations  to  yield  the  simulated  MC-MCD  spectrum  for  the  single  substitutional 
Na/Ar  system.  The  MC-MCD  simulations  of  Na/Ar  give  the  characteristic  triplet  MCD 
spectrum  with  the  correct  Boltzmann  temperature  dependence.  The  simulated 
MC-MCD  spectrum  correctly  inverts  when  the  direction  of  the  applied  magnetic  field  is 
reversed.  Addition  of  the  LCP  and  RCP  absorbances  gives  rise  to  a  characteristic  2S  — * 
2P  triplet  absorption  feature.  ©  2005  Wiley  Periodicals,  Inc.  Int  J  Quantum  Chem  103: 

854-865,  2005 


Introduction 

Alkali  metal/rare  gas  solids  (M/Rg  solids)  con¬ 
stitute  an  important  class  of  model  systems  of 
interest  to  the  U.  S.  Air  Force  in  simulating  the 
properties  of  enhanced-performance  cryogenic 
rocket  propellants  [1-4].  In  particular,  the  alkali 
metal/rare  gas  (M/Rg)  systems  share  many  phys¬ 
ical,  chemical,  and  spectroscopic  characteristics 
with  M/hydrogen  systems,  where  in  a  potential 
rocket  propulsion  application,  M  is  a  high-energy 
density  material  (HEDM),  e.g.,  a  low  mass  alkali 
metal  or  Group  13  or  14  element  [5,  6].  Boatz  and 
Fajardo  [7]  recently  modeled  2S  — >  2P  electronic 
absorption  spectra  in  M/Rg  systems  by  using  quan¬ 
tum  mechanical  first-order  degenerate  perturbation 
theory  in  conjunction  with  a  classical  Monte  Carlo 
(MC)  method  to  account  statistically  for  the  effects 
of  the  Rg  perturbations  on  the  2S  and  2P  terms  of 
the  M/Rg  systems.  Their  simulations,  while  not 
completely  quantitative,  do  recover  many  of  the 
key  attributes  of  the  electronic  absorption  spectra  of 
these  systems,  including  the  characteristic  three- 
peaked  or  "triplet"  absorption  band  shape.  Law¬ 
rence  and  Apkarian  [8]  take  a  very  similar  approach 
in  modeling  the  2P  terms  of  halogen  atoms  doped  in 
cryogenic  matrices  (X/Rg).  Boatz  and  Fajardo  make 
the  reasonable  assumption  that  spin-orbit  coupling 
represents  a  small  perturbation  of  the  2P  term  com¬ 
pared  with  the  Rg  perturbation,  at  least  for  the 
lighter  alkali  metals  of  interest  to  them  as  FIEDM 
rocket  propellant  dopants,  and  on  this  basis  justify 
the  exclusion  of  2P  spin-orbit  coupling  terms  in 
their  model.  The  Boatz  and  Fajardo  first-order  per¬ 
turbation  matrix  for  the  2P  term  is,  therefore  a  3  X 
3  matrix  whose  matrix  elements  are  defined  with 
respect  to  a  suitable  zeroth-order  orbital  p  basis  in 


the  angular  momentum  representation  (| n  l  rrij )  = 

| n  1  1),  | n  1  0),  | n  1  —1)).  Spin  is  not  considered. 
Lawrence  and  Apkarian  do  include  a  2P  spin-orbit 
term  for  the  halogen  in  their  treatment.  They  argue 
convincingly  that  the  large  halogen  2P  spin-orbit 
coupling  interaction  is  not  greatly  affected  by  the 
matrix  environment;  i.e.,  the  use  of  the  atomic 
spin-orbit  coupling  constant  of  the  halogen  in  the 
formalism  represents  a  good  approximation.  The 
Lawrence  and  Apkarian  2P  perturbation  matrix  is 
a  6  X  6  matrix  whose  individual  matrix  elements 
are  defined  with  respect  to  a  suitable  zeroth-order 
Cartesian  p  basis  with  spin  {I  nlesms)  =  I  nix—), 

I  7  1  1\  I  7  1  1\  I  7  1  1\  I  7  1  1\  I  7  1  1\1  22 

\nly-  ->,  |  n  lz--),\nlx  --->,  |  n  l  y  --->,|  n  l  z  -—>}. 

Magnetic  circular  dichroism  (MCD)  spectros¬ 
copy,  which  measures  the  differential  absorption  of 
left  circularly  polarized  (LCP)  light  vs.  right  circu¬ 
larly  polarized  (RCP)  light  in  a  sample  placed  in  a 
magnetic  field  aligned  parallel  to  the  propagation 
direction  of  the  light,  has  a  long  and  venerable 
history  of  providing  insights  to  the  nature  of  the 
chemical  and  physical  environment  surrounding 
the  MCD-active  chromophore  [9].  In  particular,  the 
MCD-active  2S  — >  2P  electronic  transitions  of  alkali 
halide  F  centers  [10-14]  and  M/Rg  systems  [15-18] 
have  attracted  considerable  experimental  and  the¬ 
oretical  attention  over  the  years.  The  usual  ap¬ 
proach  is  to  extract  parameters  such  as  the  y- value, 
the  spin-orbit  coupling  constant,  and  linear  vibra¬ 
tional  coupling  constants  from  a  moment  analysis 
of  the  experimental  MCD  spectrum  (or  spectra)  [9, 
11].  Only  a  few  full  MCD  spectral  simulations  exist 
in  the  literature.  These  simulations  of  2S  — >  2P  MCD 
spectra  rely  on  dynamic  Jahn-Teller  models  in 
which  an  a  priori  choice  is  made  as  to  which  spe¬ 
cific  interaction  vibrational  modes  will  act  to  lift  the 
degeneracy  of  the  2P  term  [19,  20].  This  simulation 
approach  is  of  limited  value  in  modeling  the  MCD 
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spectra  of  M  atoms  trapped  in  novel  and  heretofore 
unobserved  sites  in  Rg  matrices  for  which  specific 
site  geometries  and  lattice  mode  types  are  as  yet 
unknown.  A  MC  approach  for  handling  Rg  pertur¬ 
bations  in  MCD  simulations  of  M/Rg  systems  nat¬ 
urally  suggests  itself.  A  full  development  of  the 
theory  of  Monte  Carlo-Magnetic  Circular  Dichro- 
ism  (MC-MCD)  simulations  of  2S  — >  2P  MCD  spec¬ 
tra  is  presented  in  the  sections  that  follow.  The 
MC-MCD  simulation  method  is  a  first-order  de¬ 
generate  perturbation  method  that  straightfor¬ 
wardly  extends  the  MC  absorption-simulation 
methods  of  Boatz  and  Fajardo  [7]  and  Lawrence 
and  Apkarian  [8]  by  (i)  including  both  spin-orbit 
and  Zeeman  perturbations  in  the  Hamiltonian,  and 
(ii)  separately  treating  the  2S  — >  2P  transitions  asso¬ 
ciated  with  LCP  and  RCP  light.  The  MC-MCD  per¬ 
turbation  matrix  for  the  2P  term  is  a  6  X  6  matrix 
whose  individual  matrix  elements  are  expressed  in 
terms  of  a  suitable  uncoupled  angular  momentum  p 
basis  set  with  spin 

{| n  l  m,  s  ms)  =  \n  11  \  \n  1  0 ||),  \n  1 

M  |«  1  0 1  - 1>,  |n  1  -1 1  — 1>}. 

Theory 

As  noted  above,  the  theory  of  MC-MCD  spec¬ 
troscopy  of  2S  — >  2P  excitations  in  M/Rg  systems 
represents  an  extension  of  the  formalism  for  MC- 
absorption  spectroscopy  developed  by  Boatz  and 
Fajardo  [7]  (for  M/Rg  systems)  and  by  Lawrence 
and  Apkarian  [8]  (for  X/Rg  systems).  Key  aspects 
of  both  the  Boatz  and  Fajardo  and  Lawrence  and 
Apkarian  notational  schemes  are  adopted  in  the 
following  formal  development. 

TOTAL  HAMILTONIAN  AND  ZEROTII-ORDER 

EIGENVECTORS/EIGENVALUES 

The  total  Hamiltonian  for  an  alkali  metal  M  atom 
imbedded  in  a  cryogenic  rare  gas  matrix  composed 
of  N  Rg  atoms  in  the  presence  of  an  applied  exter¬ 
nal  magnetic  field  is  given  by 

H  =  Hm( r)  +  Hint( r;  R2,  R2, . .  . ,  RN),  (1) 

where 

Hm(  r)  =  Tm(  r)  +  VM(r)  (2) 


is  the  one-electron  Hamiltonian  for  the  optically 
active  electron  on  M  at  position  r,  composed  of  a 
kinetic  energy  component  Tm(t)  and  an  electrostatic 
potential  energy  component  VM(r).  The  potential 
VM(r)  is  the  potential  experienced  by  the  optically 
active  electron  in  the  free  gas-phase  M  atom  in  the 
absence  of  spin-orbit  coupling  or  an  external  mag¬ 
netic  field. 

To  describe  completely  the  2S  ground  term  of  M, 
two  zeroth-order  eigenvectors  !|«0  0^/«s);  ms  = 
1/2,  —1/2}  are  needed;  for  the  2P  excited  term,  six 
zeroth-order  eigenvectors  {| n  l  nijS  ms);  ml  .  =  1,  0, 
—  1,  ms  =  1/2,  —1/2}  are  required.  An  uncoupled 
angular  momentum  representation  {|  n  l  ml  s  ms)}  for 
the  zeroth-order  eigenvectors  is  chosen,  without 
loss  of  generality,  because  of  the  transparent  way  in 
which  MCD  selection  rules  may  be  expressed  in 
terms  of  this  representation  (vide  infra).  However, 
it  is  equally  valid  to  use  zeroth-order  eigenvector 
sets  expressed  in  the  coupled  angular  momentum 
representation  j|  n  l  s  j  in)  or  the  Cartesian  represen¬ 
tation  j| n  l  e  s  ms)}  ( e  =  x,  y,  z  for  1  =  1),  since  the 
uncoupled,  coupled,  and  Cartesian  representations 
are  mutually  related  by  unitary  transformations 
[21].  In  this  treatment,  the  interaction  Hamiltonian 
is  composed  of  four  terms,  each  of  which  depends 
on  the  nuclear  position  vectors  {R/c}  of  all  of  the  N 
Rg  atoms  in  the  matrix;  a  metal-rare  gas  interaction 
term,  a  rare  gas-rare  gas  interaction  term,  a  spin- 
orbit  term,  and  a  Zeeman  term, 

Ri,  R2,  ■  ■  ■ ,  Rn)  =  ^M-Rg(rl  R|/  IC/  •  •  •  /  Rv) 

+  KRg-Rg(r;  Ri,  R2,  ■  ■  ■  /  Rn)  +  Hso(r;  Rlr  R2, . .  . ,  Rn) 

+  ^Zeeman (h  R/  R/  •  ■  •  /  R/)-  (3) 

The  M  nucleus  is  taken  to  be  the  origin  of  the 
coordinate  system  with  respect  to  which  r  and  the 
{Rfc}  are  defined.  In  this  treatment,  it  is  assumed, 
following  Lawrence  and  Apkarian,  that  M-Rg  in¬ 
teractions  do  not  affect  Rg-Rg  interactions  [8];  no 
modifications  have  been  made  to  the  Rg-Rg  poten¬ 
tials  to  account  for  perturbations  arising  from  the 
presence  of  the  M  atom. 


Rg-Rg  AND  M/Rg  POTENTIALS 

The  Rg-Rg  interaction  potential  in  Eq.  (3)  can  be 
approximated  as  a  simple  sum  over  the  appropriate 
X  symmetry  Rg-Rg  diatomic  potentials 
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^Rg-Rg(Rl/  R-2/  •  ■  ■  /  Rn)  —  E  E  VRg_Rg(|R,„  Rj-]) 

m=l,N  m<k 

=  ^Rg-Rg  (4) 

following  the  notation  of  Boatz  and  Fajardo  [7].  It 
should  be  noted  that  URg__Rg  has  the  same  value  for 
both  the  2S  and  2P  terms  at  this  level  of  approxima¬ 
tion. 

The  M-Rg  potential,  which  in  general  exhibits  a 
complicated  dependence  on  the  Rg  positions  (Rfc),  is 
expressed  in  terms  of  a  Legendre  polynomial  ex¬ 
pansion  as 

V M-Rgh;  R1;  R2, . . . ,  Rn) 

=  E  E  Y,,(r,  R/c)  Pl(6,  (Jr,  0k,  $*),  (5) 

k=l,N  L= 0,™ 

where  r  =  ( r ,  6,  <b)  is  the  optically  active  electron 
position  vector,  and  Rfc  =  (Rk,  6k,  <f>k)  is  the  position 
vector  for  Rg  nucleus  k.  VL(r,  Rfc)  is  the  radial  po¬ 
tential  function  of  order  L,  and  PL(0,  cf>;  9k,  (f>k)  is  the 
Legendre  polynomial  of  order  L 

PL(0,  (l>;  ek,  (])k)  =  4tt(2L  +  1)~2 

x  E  rLM(e,  cf,)Y*LM(ek,  <M  (6) 

M——L,L 

and  the  YLM  are  spherical  harmonic  functions  ex¬ 
pressed  in  the  standard  Condon  and  Shortley  phase 
convention  [22].  Equation  (6)  is  an  expression  of 
one  of  the  possible  forms  of  the  spherical  harmonic 
addition  theorem  [21]. 

SPIN-ORBIT  COUPLING  OPERATOR 

The  spin-orbit  coupling  operator  for  the  optically 
active  electron  may  be  written,  using  the  standard 
defining  equation  for  the  spin-orbit  interaction  [23], 
in  terms  of  the  potential  gradient  of  the  free  M  and 
M-Rg  interaction  potentials 

Tfso(r;  Ri,  R2,  •  ■  ■ ,  Rn)  =  Hsom(t) 

+  Ytso,M-Rg(r/  Ri,  R2/  •  •  ■  /  rn)  =  {(2ffz2c2r)  1V[l/M(r) 

+  L M-Rgh/  Ri/  R2/  •  •  ■  /  Rn)]  X  p}  •  S 
=  [fiwW  +  £M-RgW]L  •  S  +  [(2mzc2r)_1 

Xv|e  E  Y,(r,  Rd  P,(0,  c/j;  9k,  X  p]  •  S  (7) 

I  k=l,N  L= l,oo 


where  p  and  S  are,  respectively,  the  linear  momen¬ 
tum  and  spin  operators  for  the  optically  active  elec¬ 
tron,  L  =  r  X  p  is  the  angular  momentum  operator 
for  the  optically  active  electron,  and 

£M(r)  =  (2m2c2r)_1d[  VM(r)]/ dr  (8) 

is  the  standard  spin-orbit  coupling  operator  for  the 
free  M  atom.  The  M-Rg  interaction  potential  of  Eq. 
(5)  can  be  broken  up  into  spherical  and  nonspheri- 
cal  components 

l^M-RgCc  Ri/  Rr/  •  •  •  /  Rn)  =  E  V0(r,  Rjt) 

k—l,N 

+  E  E  VL(r,  IWPl(0,0;  Bh  4).  (9) 

k=l,N  L— 1,0° 

In  the  form  of  the  theory  currently  in  use  for 
MC-MCD  simulations,  it  is  assumed  that  the  gra¬ 
dient  expression  in  Eq.  (7)  is  small  and  can  be 
neglected;  i.e.,  no  L  >  0  terms  are  assumed  to 
contribute  to  the  spin-orbit  coupling  interaction. 
The  total  spin-orbit  coupling  operator  of  the  M/Rg 
system  in  the  spherical  approximation  is  thus  given 
by 

Yfso(r)  =  HSOM(t)  +  HSOiM_Rg(r) 

=  [Ur)  Y  &,_Rg(r)]L  '  S  =  &T,  L  •  S),  (10) 

where  the  dependence  on  the  Rg  position  vectors 
[Rfc],  while  present,  is  not  explicitly  shown  in  the 
interest  of  preserving  a  compact  notation  in  the 
final  form  of  the  equations  used. 

No  spin-orbit  coupling  is  possible  when  the  op¬ 
tically  active  electron  is  in  the  2S  ground  state  man¬ 
ifold,  where  l  =  0.  In  the  2P  excited  state  manifold, 
where  Z  =  1,  the  spin-orbit  coupling  constant  is 
expressed  formally  as  the  sum  of  the  radial  expec¬ 
tation  values  of  the  spherically  symmetric  M  and 
M-Rg  spin-orbit  coupling  operators 

Cni  =  (n  l\U)\n  l >  =  <■ n  l\%M(r)\n  l ) 

+  (n  Z|^M-Rg(r)ln  Z)  =  £m  +  £m-rg-  (11) 

If  the  radial  eigenvector  | n  l)  is  the  free  atom  P  state 
atomic  radial  eigenvector  for  M,  or  if  | n  l )  closely 
approximates  this  eigenvector,  then  is  the  free 
alkali  atom  spin-orbit  coupling  constant  and  fM_RG 
is  the  spherical  perturbation  to  the  spin-orbit  cou¬ 
pling  constant  induced  in  M  by  the  Rg  matrix  en¬ 
vironment.  Equation  (11)  provides  an  insightful 
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context  for  discussing  the  possibilities  of  (1)  spin- 
orbit  coupling  constant  sign  reversal  or  (2)  reduc¬ 
tion  of  the  magnitude  of  a  spin-orbit  coupling  con¬ 
stant  (partial  quenching),  in  a  2P  M-Rg  system.  For 
any  given  alkali  metal  M,  the  free  atom  2P  spin- 
orbit  coupling  constant  is  positive;  i.e.,  £M  >  0. 
Thus,  for  a  2P  spin-orbit  coupling  constant  £n/  to  be 
negative  in  a  M-Rg  system,  the  M-Rg  spin  orbit 
contribution  must  be  negative  and  larger  in  magni¬ 
tude  than  the  free  M  atom  spin-orbit  coupling  con¬ 
stant;  i.e.,  £m_rg  <  0  and  |£M_RG|  >  £M.  The  condi¬ 
tions  for  the  partial  quenching  case,  where  £nl  is  still 
positive  but  smaller  than  £M,  are:  £M-rg  <  0  and 
|£m-rgI  <  £m-  in  the  actual  MC-MCD  simulations, 
u  is  introduced  as  an  adjustable  parameter.  The 
two  cases  may  be  explored  in  the  MC-MCD  simu¬ 
lations  by  making  appropriate  choices  for  £n,  with 
respect  to  the  known  positive  £m  of  the  alkali  metal 
M  atom  of  interest. 

ZEEMAN  PERTURBATION 

The  Zeeman  perturbation  operator  is  given  by 

i  i  /cerium  T/  R-l/  R-2/  *  *  •  /  R \  i  —  AB^localCL  +  2S),  (12) 

where  pB  is  the  Bohr  magneton,  L  and  S  are  the 
orbital  and  spin  angular  momentum  operators  for 
the  optically  active  electron,  and  Blocal  is  the  mag¬ 
netic  field  locally  experienced  by  the  optically  ac¬ 
tive  electron  as  a  result  of  the  application  of  an 
external  magnetic  field  B  to  the  M/Rg  system.  Both 
the  magnitude  and  the  direction  of  the  local  mag¬ 
netic  field  vector  are  influenced  by  the  positions  of 
the  Rg  atoms  and  by  the  nature  of  the  electronic 
state  of  the  optically  active  electron  (i.e.,  2S  or  2P). 
The  relationship  between  the  external  magnetic 
field  and  the  local  magnetic  field  thus  may  be  ex¬ 
pressed  as 

Blocal  B.v] ! (  R  ; ,  R2,  •  •  •  ,  Rn) 

=  [1  -  aJR„  R2 . Rn)]B,  (13) 

where  <r„;(Ri,  R2,  ■ . . ,  R:V)  is  a  3  X  3  shielding  tensor 
whose  individual  components  depend  on  both  the 
nature  of  the  electronic  state  being  shielded  and  the 
positions  of  the  Rg  atoms  in  the  matrix.  In  actual 
MC-MCD  simulations,  it  is  assumed  that  the  local 
magnetic  field  is  aligned  along  the  laboratory  z-axis 

R  local  =  f>localez/  (14) 


which  gives 

^Zeeman  =  AbB(1  “  n„,)(Lz  +  2SZ)  (15) 

Different  values  of  Blocal  can  be  chosen,  if  desired,  to 
model  the  effects  of  slightly  different  2S  and  2P 
shielding  constants  an0  and  anl.  Changes  in  the 
direction  of  Blocal  away  from  the  z-axis  arising  from 
the  tensor  character  of  crnl  can  be  modeled  by  pre¬ 
serving  the  laboratory  z-axis  orientation  of  the  Zee- 
man  perturbation,  which  is  very  desirable  from  the 
standpoint  of  matrix  element  evaluation,  and  rotat¬ 
ing  the  Rg  position  vector  set  through  an  arbitrary 
unitary  rotation  transformation  U{Rk)  =  |R(.(.  In  any 
event,  it  is  expected  that  differences  between  B  and 
Biocai  will  be  comparatively  small  for  closed  shell 
Rg  atoms  in  their  ground  S  electronic  states.  Cer¬ 
tainly  within  the  chosen  MC-MCD  simulation  re¬ 
strictions  mentioned  earlier  of  neglecting  polariza¬ 
tion  effects  and  the  higher  excited  states  of  the  Rg 
atoms,  B  will  not  differ  much  from  Blocal. 


DIAGONALIZATION  OF  TIIE  FULL 
HAMILTONIAN  MATRIX  IN  THE  2S  AND  2P 
MANIFOLDS 

An  implementation  of  degenerate  first-order  per¬ 
turbation  theory  is  required  to  arrive  at  approxi¬ 
mate  2S  and  2P  eigenvalues  and  eigenvectors  of  the 
full  M/Rg  Flamiltonian  given  in  Eq.  (1).  These  eig¬ 
envalues  and  eigenvectors  contain  the  effects  of  the 
various  terms  in  the  interaction  Flamiltonian  Eq. 
(3).  First-order  perturbation  theory  gives  eigenval¬ 
ues  that  are  correct  through  first  order.  The  associ¬ 
ated  eigenvectors  are  those  particular  zeroth-order 
eigenvectors,  expressed  as  orthonormal  linear  com¬ 
binations  of  the  zeroth-order  eigenvector  set 
{\n  I  m,  s  ms)}, 

^i=\nli)=  2  X  cfmism’\n  1  m,  s  ms), 

mi—  —  1,1  ms  =  - 1/2,1/ 2 

(16) 

(n  1  i'\n  1  i )  =  8(P,  i),  (17) 

that  diagonalize  the  2  X  2  (2S  manifold)  or  6  X  6  (2P 
manifold)  matrix  formed  by  sandwiching  the  full 
Hamiltonian,  Eq.  (1),  within  the  same  uncoupled 
zeroth-order  eigenvector  set  {| n  I  m,s  ms)};  i.e.. 
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(n'l'm'is'm's\H\n  l  m,s  ms ) 

=  (n1  V m[s' m's\HM\n  l  mt  s  ms) 

+  (n'rm'is'm's\Hint\n  1  mt  s  ms ).  (18) 

Both  the  first-order  eigenvalues  and  the  associated 
zeroth-order  eigenvectors  are  needed  for  MC-MCD 
simulations.  Since  the  first  matrix  element  on  the 
right  side  of  Eq.  (18)  is  the  zeroth-order  matrix 
element  of  the  free  M  Hamiltonian,  it  is  automati¬ 
cally  diagonal  in  either  the  original  zeroth-order 
basis  [\nlmjsmJ)} 

( n'rm'is'm's\HM\n  l  mts  m$) 

=  e„,8(n' ,  n)8(l' ,  l)8{m\,  m,)8{s' ,  s)8(m's,  ms) 

(19) 

or  the  new  zeroth-order  basis  j| n  li)} 

( n  1  i'\HM\n  1  i)  =  e„,8(n',  n)8(l',  l)8(i',  i ).  (20) 

Within  the  2S  manifold,  the  interaction  matrix  on 
the  right  side  of  Eq.  (18)  also  is  diagonal;  thus 

dy(2S)  =  \n  0  i)  =  2  Cf0m‘\n  0  0  \  ms), 

ms=—l/ 2,1/  2 

(21) 

which  reduces  to 

dh(2S)  =  \n  0  0  |f>,  (22a) 

¥2(2S)  =  \n  0  0  |  (22b) 

The  i  index  is  mapped  onto  the  ms  index  one-to- 
one  in  this  specific  case.  The  final  expressions  for 
the  2S  eigenvalues,  correct  through  first  order,  are 
obtained  by  collecting  together  the  matrix  elements 
previously  derived  for  the  Rg-Rg,  M-Rg,  and  Zee- 
man  perturbations: 

E;(2S)  =  e„0(2S)  +  En0i(2S)  =  e„0(2S)  +  ltRg_Rg(2S) 

+  UM_Rg(~S)  +  LtZeeman(2S) 

=  eiio(2S)  +22  ^Rg-Rg/XS ( I R/ji  —  R*|) 

m—l,N  m<k 

+  22  ^M-Rg,Xs(R)  +  AflBloca l(2 ws) .  (23) 

m—l,N  m<k 


In  Eq.  (23),  the  eigenvalue  convention  of  Boatz 
and  Fajardo  [7]  is  adopted:  £,  is  an  eigenvalue  of  a 
specific  2S  state  correct  through  first  order,  e„0  is 
the  zeroth  order  2S  eigenvalue  for  the  optically 
active  electron  of  the  free  M  atom,  and  En0i  is  the 
first-order  eigenvalue  correction  term  arising  from 
the  interaction  Hamiltonian. 

To  obtain  2P  final  state  eigenvalues  correct 
through  first  order 

Ej{2P )  =  snl(2P )  +  Enl)(2P),  (24) 

the  first-order  eigenvalue  set  (Enlj(2P);f  =  1, 6)  must 
be  calculated  by  diagonalizing  the  6  X  6  2P  inter¬ 
action  matrix  in  Eq.  (18).  The  individual  matrix 
elements  of  this  interaction  matrix,  shown  in  Table 
I,  are  obtained  by  collecting  together  the  matrix 
elements  previously  derived  for  the  2P  Rg-Rg, 
M-Rg,  spin-orbit,  and  Zeeman  perturbations.  It 
should  be  noted  that  the  l/Rg_Rg  and  the  (l/0)M_Rg 
perturbations  are  diagonal  perturbations  in  this  6  X 
6  matrix  and  can  be  separated  from  the  6X6  matrix 
eigenvalue  problem  as  illustrated  by  Boatz  and  Fa¬ 
jardo  [7]  in  their  Eqs.  (6)  and  (8).  The  associated 
final  state  2P  eigenvectors  that  diagonalize  this  ma¬ 
trix,  are 

d'/P)  =  \nlf) 

=  2  2  C;lmi(1/2)>  1  mx\ma),  (25) 

mi=  — 1,1  ms=  — 1/  2,1/  2 

where  /  =  1,6  and 

(n  1  f\n  1  /)  =  8(f, /).  (26) 

2S  ->  2P  TRANSITION  ENERGIES  AND 

TRANSITION  MOMENT  INTEGRALS  FOR 

M/Rg  MCD  SPECTROSCOPY 

The  combined  effects  of  the  Rg-Rg,  M-Rg,  spin- 
orbit,  and  Zeeman  perturbations  on  the  optically 
active  electron  of  a  M/Rg  system  act  to  split  the 
ground  2S  manifold  into  two  energetically  distinct 
states  {d',;  i  =  1,  2}  and  the  excited  2P  manifold  into 
six  energetically  distinct  states  { /  =  1,  6}.  The 
specific  characteristics  of  the  2S  and  2P  eigenvalues 
and  their  associated  eigenvectors  are  dependent 
upon  the  particular  Rg  configuration  { R, ,  R2, .  .  . , 
Rn)  within  which  the  problem  is  solved.  Thus,  for  a 
given  Rg  configuration  {R1,  R2, .  .  . ,  RN)  and  mag¬ 
nitude  of  the  magnetic  field,  up  to  12  energetically 
and  optically  distinct  transitions  may  be  identified 
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TABLE  I  _ 

Matrix  elements  of  the  6x6  2P  interaction  matrix  in  the  uncoupled  angular  momentum  representation  |n/m^ms). 

|n1iy>  |n10y>  Ini  —  1  |n11  \  ~  |n10 
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\V0)k  =  j  [VBX(Rk)  +  2VAn{Rk)]  and  {V2)k  =  5[ VBS{Rk )  -  VAn(Rk)].  The  final  matrix  elements  are  constructed  by  summing  over  “k”,  which  indexes  the  rare  gas  atoms. 
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within  the  manifold-to-manifold  2S  — >  2P  transition 
of  the  M/Rg  system.  The  transition  energies  of 
these  12  2S  ~^2P  transitions  are  given  as 

hvfi  =  Ef(2P)-EieS),  (27) 

where  i  =  1,  2  and  /  =  1,  6.  The  Rg-Rg  interaction 
energy  terms  cancel  in  Eq.  (27),  since  URg_Rg(2P)  = 

^Rg-Rg(  S)  =  Ft Rg— Rg- 

Electronic  transitions  involving  the  optically  ac¬ 
tive  electron  of  a  M/Rg  system  are  well  described 
as  electric  dipole  transitions.  In  MCD  spectroscopy, 
each  of  the  12  transitions  may  be  induced  to  a 
greater  or  lesser  extent  by  either  LCP  or  RCP  light 
interacting  with  the  M/Rg  system  in  the  electric 
dipole  limit  and  tuned  to  the  particular  transition 
energy  hvc.  Thus,  both  LCP  and  RCP  electric  dipole 
transition  moment  integrals  must  be  computed  for 
each  of  the  12  possible  combinations  of  initial  and 
final  states;  i.e.,  24  transition  moment  integrals,  12 
for  LCP  light  and  12  for  RCP  light,  are  needed  to 
describe  the  MCD  spectrum  of  a  2S  — >  2P  transition 
in  M/Rg  system  in  a  particular  Rg  configuration 
{Ri,  R2/  ■  •  ■  /  Rn1;  i-ev 

<^/2P)|0LCp|^(2S)>  (28a) 

(^/2P)|Orcp|^(2S)),  (28b) 

where  i  =  1,2  and  /  =  1,  6.  Olcp  and  Orcp  are  the 
electric  dipole  transition  operators  for  LCP  and 
RCP  light  in  Eqs.  (28a)  and  (28b).  The  electric  dipole 
moment  operator  m  of  the  M/Rg  system  is,  to  a 
very  good  initial  approximation,  determined  by  the 
position  r  =  (x,  y,  z)  and  charge  —  e  of  the  optically 
active  electron.  The  Rg  atoms  are  electrically  neu¬ 
tral,  and  their  contributions  to  the  dipole  moment 
can  be  neglected  if  polarization  effects  are  assumed 
to  be  small.  The  M+  core  is  positioned  at  the  origin 
of  the  chosen  coordinate  system  and  thus  does  not 
contribute  to  the  expression  for  the  dipole  moment 
at  this  level  of  approximation,  i.e., 

m  =  X  cji'i'i  =  +  eO  =  -e(xex  +  yey  +  zez) 

i 

+  e(0ex  +  0ev  +  0e:)  =  —  e(xex  +  yey  +  zez).  (29) 

However,  the  optically  active  electron  at  position 
r  will  induce  a  slight  polarization  in  the  M+  core 
and  the  Rg  matrix.  This  polarization  can  be  repre¬ 
sented  as  a  small  negative  electrostatic  image 


charge  —q’  at  —  r.  Thus,  polarization  effectively  re¬ 
duces  the  charge  of  the  optically  active  electron  to 

— e'  =  —e  +  q' ,  (30) 

where  \—e'\  <  |— e\,  and  the  dipole  moment  operator 
of  the  M/Rg  system  becomes 

m  =  — e'r  =  —e'(xex  +  yey  +  zez).  (31) 

While  it  can  be  argued  that  the  image  charge  —q' 
(and  hence  —e')  is  not,  strictly  speaking,  a  constant, 
but  will  exhibit  both  a  radial  dependence  (i.e.,  dif¬ 
ferent  q'  for  small  r  vs.  large  r)  and  angular  depen¬ 
dence  arising  from  the  specific  positions  of  the  Rg 
atoms  (R],  R2, . . . ,  R;V)  in  the  M/Rg  system,  these 
effects  are  neglected  at  this  level  of  the  theory. 
Within  the  Piepho  and  Schatz  unit  vector  conven¬ 
tions  for  LCP  and  RCP  light  [9], 


'n'Lep  =  2  1/2(ex  -  zev) 

(32a) 

ttrcp  =  2“1/2(ex  +  zey). 

(32b) 

the  LCP  and  RCP  electric  dipole  moment  operators 
for  the  M/Rg  system  become 

Olcp  =  m  ■  tt*Cp  =  ~e'Tl,\x  +  iy) 

=  e'(477/3)1/2Yn(0,  4>)r,  (33a) 

Orcp  =  m  •  tt|cp  =  -e'Tx/\x  -  iy) 

=  — e'(4'7r/3)1/2Y1_1(0,  <«r.  (33b) 

Since  each  of  the  ^P,(2S)  and  'I,/2P)  states  is  ex¬ 
pressed  in  terms  of  the  uncoupled  zeroth-order  an¬ 
gular  momentum  eigenvector  set  j| nlmpm^}  [see 
Eqs.  (16)  and  (25)],  the  specific  form  of  the  transition 
moment  integrals  given  in  Eqs.  (33a)  and  (33b)  may 
be  the  determined  by  evaluating  integrals 

(n  Vm  Is '  m  ( |  O LCP|  n hiy  sms) 

=  ( nl'm ! \  ms'|OLCP|n00  \  ms)  =f8(ml,  1  )S(ms',  ms), 

(34a) 

(n  l’m\  s'm's\0RCp\nlmisms) 

1 

=  {nl'm'i  ^  /«s'|ORCP]«00  \  ms) 

=  —T)8(m's,  ms),  (34b) 
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where 

/lcp  =/=  (e'[4'7r/3]1/2)<nl|r|nO)<Y11|Y11|Y00) 

=  (e'[3Y1/2)(nl\r\nO)  (35a) 

/rcp  =  -/=  (-e'^ir/Sp/^^lHnOXYnlY^lYoo) 

=  (e'[3]~1/2)(nl|r|«0)  (35b) 

Equations  (34a)  and  (35a)  express  the  fact  that  an 
LCP  electric  dipole  transition  acts  on  the  orbital 
angular  momentum  component  of  the  uncoupled 
wave  function,  increasing  both  /  and  mt  by  one  unit 
while  leaving  s  and  ms  unchanged.  The  correspond¬ 
ing  RCP  electric  dipole  transition  described  in  Eqs. 
(34b)  and  (35b)  also  acts  on  the  orbital  angular 
momentum  in  the  uncoupled  wave  function,  in¬ 
creasing  l  by  one  unit  and  decreasing  ml  by  one 
unit.  In  Eqs.  (35a)  and  (35b), /is  the  magnitude  of 
the  evaluated  transition  moment  integral.  Note  that 
the  value  of/ is  the  same  for  all  allowed  transitions. 
By  combining  the  transition  moment  integral  ex¬ 
pressions  in  Eqs.  (28),  (34),  and  (35)  with  the  appro¬ 
priate  expressions  for  T',(2S)  and  d/(2P)  in  the  un¬ 
coupled  basis  as  given  in  Eqs.  (16),  (17),  and  (25), 
final  expressions  for  the  contribution  to  the  MCD 
spectrum  by  the  transition  i  — »/ at  energy  h  ty  can  be 
written: 

A  A'(vfi)  =  ALCP(iy)  —  ARCP(vfi) 

=  a{|(^/2P)|0LCP|^i(2S))|2  -  |<^/2P)|Orcp|^,(2S))|2} 

_  qC 2j|£Ml,m/=l,ms|2  |^nl/m/=-l/ms|2| 

X  [1  +  (e~2^B/kT  -  1  )8(ms,  1/2)],  (36) 

where  a  is  a  numerical  constant  that  relates  the 
squared  transition  moment  integrals  to  the  absor¬ 
bances  and  the  initial  state  index  i  is  associated  with 
ms  =  —1/2  or  1/2  in  the  2S  manifold.  The  MCD 
spectral  amplitude  at  energy  hvc  for  the  Rg  config¬ 
uration  {Rj,  R2, .  . . ,  Rn}  is  seen  to  be  related  to 
differences  of  squares  of  selected  2P  eigenvector 
coefficients  computed  for  that  configuration.  The 
term  in  square  brackets  on  the  right  side  of  Eq.  (36) 
is  a  Boltzmann  factor  for  the  2S  manifold.  The  Zee- 
man  perturbation  splits  the  energies  of  the  ms  = 
—  1/2  and  1/2  states  of  the  2S  manifold  by  a  factor 
of  2/rBB.  The  ms  =  —1/2  and  1/2  states  will,  there¬ 
fore,  have  different  thermal  populations.  This  must 
be  accounted  for  in  computing  the  various  contri¬ 
butions  to  the  MCD  spectrum. 


Monte  Carlo  Simulation  of  the  MCD 
Spectrum 

The  full  MCD  spectrum  for  a  2S  — >  2P  transition 
in  a  M/Rg  system  in  a  specific  Rg  configuration  { R , , 
R2, .  . . ,  R;V1  is,  therefore,  a  "stick  MCD  spectrum" 
consisting  of  12  lines,  appropriately  placed  on  the  x 
(energy)  axis  at  transition  energies  hvfi/  whose  mag¬ 
nitudes  and  directions  (+  or  — )  on  the  y  (A A’)  axis 
are  determined  by  computing  Eq.  (36)  for  all  pos¬ 
sible  choices  of  i  and  /.  To  simulate  a  real  MCD 
spectrum,  many  stick  MCD  spectra,  each  arising 
from  a  different  Rg  configuration,  are  averaged  to¬ 
gether  using  a  simple  binning  (histogram)  tech¬ 
nique  based  on  the  original  Metropolis  et  al.  [24] 
algorithm  exactly  as  implemented  by  Boatz  and 
Fajardo  [7]  for  the  case  of  electronic  absorption 
spectroscopy.  A  MC-MCD  simulation  begins  by 
choosing  an  initial  Rg  configuration  { R , ,  R2, .  .  . , 
R;\'l initial/  typically,  representing  an  idealized  trap¬ 
ping  site  (e.g.,  a  single  substitutional  site  with  Oh 
symmetry).  The  MC  energy  optimization  scheme 
operates  in  the  2S  ground  state  manifold.  At  this 
stage,  the  very  small  Zeeman  perturbation  of  the  2S 
manifold  is  neglected.  MCD  stick  spectra  arising 
from  more  favorable  configurations  are  weighted 
more  heavily  in  the  averaging  process  than  are 
those  arising  from  the  less  favorable  configurations. 
In  actual  practice,  the  LCP  and  RCP  contributions 
to  the  MCD  spectrum  appearing  in  Eq.  (36)  are 
accumulated  and  stored  separately.  This  allows  the 
electronic  absorption  spectrum  to  be  recovered 
from  the  MC-MCD  simulation  as  [9] 

A’  =  A  =  \  [A  LCp  +  ArCP].  (37) 

Equation  (37)  expresses  the  fact  that  the  electronic 
absorption  spectra  of  M/Rg  systems  are  essentially 
unaffected  by  the  application  of  an  external  mag¬ 
netic  field. 

ORIENTATIONAL  AVERAGING  IN  MC-MCD 

SIMULATIONS 

In  an  MCD  experiment,  the  propagation  direc¬ 
tion  of  the  light,  which  is  always  parallel  to  the  B 
vector,  rigorously  defines  the  laboratory  frame  z- 
axis.  The  laboratory  frame  x-  and  y-axes  are  set 
parallel  to  the  sapphire  deposition  window  and 
perpendicular  to  B  with  the  origin  at  the  M  nucleus. 
All  Rg  positions,  eigenvectors,  and  selection  rules 
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are  described  with  respect  to  this  laboratory  frame. 
To  model  the  effects  of  arbitrary  Rg  lattice  orienta¬ 
tion  relative  to  the  laboratory  frame,  the  6  X  6  2P 
matrix  listed  in  Table  I  is  separated  into  M/Rg  and 
spin-orbit  +  Zeeman  matrices: 

[M/Rg  +  SO  +  Zeeman] 

=  [M/Rg]  +  [SO  +  Zeeman],  (38) 

Arbitrary  orientations  of  the  Rg  matrix  relative 
to  the  laboratory  frame  can  be  modeled  by  rewrit¬ 
ing  Eq.  (38)  as 

[M/Rg  +  SO  +  Zeeman]  =  [U][M/Rg][U]_1 

+  [SO  +  Zeeman],  (39) 

where  U  is  an  arbitrary  randomly  generated  uni¬ 
tary  rotation  matrix  (e.g.,  a  Eulerian  rotation  ma¬ 
trix).  The  MC-MCD  simulation  process  is  carried 
out  as  before,  using  Eq.  (39)  rather  than  Eq.  (38). 

LATTICE  TEMPERATURE  AND  MAGNETIC 

TEMPERATURE  IN  MC-MCD  SIMULATIONS 

Temperature  effects  appear  in  two  distinct  places 
in  MC-MCD  simulations.  A  lattice  temperature 
must  be  specified  for  the  MC  simulation  process.  It 
may  be  argued  that  a  proper  representation  of 
quantum  mechanical  zero  point  vibrational  mo¬ 
tions  of  the  M/Rg  lattice  at  temperature  T  requires 
the  classical  MC  simulation  to  be  carried  out  at  a 
classical  lattice  temperature  T  >  T  [8].  However,  a 
temperature  for  the  Boltzmann  factor  of  the  Zee¬ 
man-split  2S  term  also  must  be  included.  This  mag¬ 
netic  temperature,  Tmag,  should  always  be  the  ac¬ 
tual  system  temperature  T  (Tmag  =  T).  The  MC- 
MCD  simulation  code  has  the  option  of  specifying 
separate  values  for  the  lattice  temperature  T'  and 
the  magnetic  temperature  Tmag.  Lawrence  and  Ap- 
karian  [8]  discuss  the  relationship  between  T'  and  T 
in  classical  MC  simulations  of  doped  Rg  systems.  In 
practice  T'  is  calibrated  to  T  by  comparing  the  Rg 
MC  equilibrium  lattice  spacing  (for  first  and  second 
nearest  neighbors)  with  the  experimentally  mea¬ 
sured  values  [25].  For  Ar  at  10  K,  the  first  nearest- 
neighbor  analysis  predicts  a  lattice  temperature  of 
31  K,  while  the  second  nearest-neighbor  analysis 
predicts  a  slightly  lower  lattice  temperature  of  27  K. 
Consequently  a  simulation  lattice  temperature  in 
the  vicinity  of  30  K  represents  the  best  modeling  of 
the  quantum  zero  point  vibrations  of  the  Ar  lattice. 
The  Na/Ar  system  at  a  magnetic  temperature. 


Simulated  Absorption  Spectra 


Simulated  MCD  Spectra 


Energy  (cm'1) 

(b) 


FIGURE  1.  Comparison  of  the  temperature  effects  in 
the  (a)  absorption,  and  (b)  MCD  spectra  for  the  case  of 
SO  =  17  cm-1,  B  =  1.0  T,  magnetic  temperature, 

Tm  =  10  K,  at  simulation  temperatures  of  10,  30,  and 
45  K. 


Tmag  =  10  K,  a  Na  2P  spin-orbit  parameter  of  17 
cm-1  (the  value  for  unperturbed  atomic  Na)  is  il¬ 
lustrated  in  Figure  l(a,  b),  in  which  are  presented 
the  simulated  absorption  and  MCD  spectra,  respec¬ 
tively,  for  three  different  choices  of  lattice  temper¬ 
ature,  T'  =  10,  30,  and  45  K. 
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FIGURE  2.  Absorption  spectra  with  (a)  normal  and  (b)  reversed  B-fields.  1-atom  vacancy,  SO  =  17  cm  1 ,  and 
Tmag  =  6,  8,  10  K.  The  magnitudes  of  the  B-fields  are  1.0  T. 


While  the  characteristic  three-peaked  "triplet" 
absorption  feature  and  the  corresponding  down- 
down- up  MCD  triplet  feature  are  preserved  for  all 
three  choices  of  T',  the  spectral  breadth  increases  as 
T'  increases,  improving  agreement  between  simu¬ 
lated  and  experimental  spectra.  The  difference  in 
peak  height  is  due  to  the  fact  that  the  total  area 
under  the  absorption  curve  is  normalized  to  1.0; 
broader  peaks  will  be  shorter. 

OTHER  IMPORTANT  PROPERTIES  OF 
MC-MCD  SIMULATIONS 

For  the  single  substitutional  site  in  the  Na/Ar 
system,  absorption  and  MC-MCD  simulations  cor¬ 
rectly  reproduce  other  observed  properties  of  ex¬ 


perimental  absorption  and  MCD  spectra  M/Rg  sys¬ 
tems.  For  example,  the  observed  insensitivity  of  the 
experimental  2S  — >  2P  M/Rg  triplet  absorption  fea¬ 
ture  to  the  magnetic  field  strength  and  Tmag  is  faith¬ 
fully  reproduced  in  the  MC-MCD  simulation,  as 
shown  in  Figure  2.  Reversal  of  the  magnetic  field 
direction  (B  =  Bez  changed  to  B  =  — Bez)  in  an  MCD 
experiment  has  the  effect  of  inverting  the  MCD 
spectrum.  If  the  normal  MCD  spectrum  follows  a 
down- down- up  pattern,  the  MCD  spectrum  with 
the  reversed  magnetic  field  will  be  the  mirror  image 
of  the  original  spectrum  with  an  up- up- down  pat¬ 
tern.  Figure  3  shows  that  the  MC-MCD  simulations 
of  Na/  Ar  correctly  give  an  inverted  spectrum  when 
B  is  replaced  with  -B.  These  successes  suggest  the 


(a) 


(b) 


FIGURE  3.  MCD  spectra  with  (a)  normal  magnetic  field  alignment,  and  (b)  reversed  magnetic  field  alignment. 
1-atom  vacancy,  SO  =  17  cmr1,  and  7mag  =  10  K.  The  magnitudes  of  the  B-fields  are  1.0  T. 
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MC-MCD  simulation  embodies  much  of  the  essen¬ 
tial  physics  of  the  Na/Ar  system  and  its  interaction 
with  external  magnetic  fields,  the  perturbing  effects 
of  the  Rg  matrix  and  electromagnetic  radiation. 
Overall,  the  MC-MCD  simulated  Na/Ar  absorp¬ 
tion  and  MCD  spectra  show  good  qualitative  agree¬ 
ment  with  experimental  Na/Ar  spectra. 

ACKNOWLEDGMENTS 

It  is  a  great  pleasure  to  acknowledge  Mario  E. 
Fajardo  for  his  helpful  discussions  on  theory  and 
for  his  insight  in  making  connections  between  the¬ 
ory  and  experiment. 


References 

1.  Fajardo,  M.  E.;  Carrick,  P.  G.;  Kenney,  J.  W.  III.  J  Chem  Phys 
1991,  94,  5812. 

2.  Fajardo,  M.  E.  J  Chem  Phys  1993,  98,  110. 

3.  Tam,  S.;  Fajardo,  M.  J  Chem  Phys  1993,  99,  854. 

4.  Corbin,  R.;  Fajardo,  M.  J  Chem  Phys  1994,  101,  2678. 

5.  Thompson,  T.  L.;  Cordonnier,  M.  E.;  Fajardo,  M.  E.;  Tam,  S. 
Chem  Phys  1994,  189,  361. 

6.  Carrick,  P.  G.  Specific  impulse  calculations  of  high  energy 
density  solid  cryogenic  rocket  propellants.  I.  Atoms  in  solid 
h2.  Technical  Report  PL-TR-93-3014;  Phillips  Laboratory:  Ed¬ 
wards  AFB,  CA,  1993. 

7.  Boatz,  J.;  Fajardo,  M.  J  Chem  Phys  1994,  101,  3472. 

8.  Lawrence,  W.  G.;  Apkarian,  V.  A.  J  Chem  Phys  1994,  101, 
1820. 


9.  Piepho,  S.  B.;  Schatz,  P.  N.  Group  Theory  in  Spectroscopy 
With  Applications  to  Magnetic  Circular  Dichroism;  Wiley: 
New  York,  1983. 

10.  Mort,  J.;  Luty,  F.  Phys  Rev  Lett  1964,  12,  45. 

11.  Henry,  C.  H.;  Snatterly,  S.  E.;  Schlickter,  C.  P.  Phys  Rev  A:  At 
Mol  Opt  Phys  1965,  137,  A583. 

12.  Brown,  F.  C.;  Mort,  J.;  Luty,  F.  Phys  Rev  A:  At  Mol  Opt  Phys 
1965,  137,  566. 

13.  Merle  d'Aubinge,  Y.;  Roussel,  A.  Phys  Rev  B:  Condensed 
Matter  and  Material  Physics  1971,  3,  1421. 

14.  Stephens,  P.  J.;  Osborne,  G.  A.  J  Chem  Phys  1972,  56,  609. 

15.  Hormes,  J.;  Schiller,  J.  Chem  Phys  1983,  74,  433. 

16.  Lund,  P.  A.;  Smith,  D.;  Jacobs,  S.  M.;  Schatz,  P.  N.  J  Phys 
Chem  1984,  88,  31. 

17.  Rose,  J.;  Smith,  D.;  Williamson,  B.  E.;  Schatz,  P.  N.;  O'Brien, 
M.  C.  M.  J  Phys  Chem  1986,  90,  2608. 

18.  Samet,  C.;  Rose,  J.  L.;  Williamson,  B.  E.;  Schatz,  P.  N.  Chem 
Phys  Lett  1987,  142,  557. 

19.  O'Brien,  M.  C.  M.  J  Phys  C:  Solid  State  Phys  1985,  18,  4963. 

20.  O'Brien,  M.  C.  M.  J  Chem  Phys  1985,  82,  3870. 

21.  Weissbluth,  M.  Atoms  and  Molecules;  Academic  Press:  San 
Diego,  CA,  1978. 

22.  Condon,  E.  U.;  Shortley,  G.  H.  The  Theory  of  Atomic  Spectra; 
Cambridge:  London,  1963. 

23.  Griffith,  J.  S.  Theory  of  Transition  Metal  Ions;  Cambridge: 
London,  1961. 

24.  Metropolis,  N.;  Rosenbluth,  A.  W.;  Rosenbluth,  M.  N.;  Teller, 
A.  H.;  Teller,  E.  J  Chem  Phys  1953,  21,  1087. 

25.  Terrill-Stolper,  H.  A.  Theory  of  the  Magnetic  Circular  Di¬ 
chroism  (MCD)  and  Electronic  Absorption  Spectra  of  Alkali 
Metals  Matrix  Isolated  in  Cryogenic  Rare  Gases.  I;  Master's 
thesis;  New  Mexico  State  University,  Las  Cruces,  NM,  1996. 


INTERNATIONAL  JOURNAL  OF  QUANTUM  CHEMISTRY 


865 


